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Abstract. We investigate the low-energy properties of a generalized quantum sine-Gordon 
model in one dimension with a self-dual symmetry. This model describes a class of quantum 
phase transitions that stems from the competition of different orders. This SU(7V) self-dual 
sine-Gordon model is shown to be equivalent to an SO(A02 conformal field theory perturbed 
by a current-current interaction, which is related to an integrable fermionic model introduced 
by Andrei and Destri. In the context of spin-chain problems, we give several realizations of 
this self-dual sine-Gordon model and discuss the universality class of the transitions. 
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Duality symmetries have been providing much insight in diverse areas of physics, 
ranging from high-energy physics to condensed matter physics or statistical mechanics. One 
of the main reasons for this is that a duality maps, in general, a theory in the strong coupling 
onto one in the weak coupling, and thus is a powerful tool for investigating strongly coupled 
regimes. In some lattice spin models, the duality transformation can be carried out explicitly, 
mapping the partition function of one theory to that of another or to the same theory if the 
theory is self-dual. The simplest well-known example is the Kramers-Wannier (KW) duality 
transformation of the two-dimensional Ising model, which even locates the critical point 
without calculating the partition function explicitly [ l j. In the context of the equivalent one- 
dimensional quantum Ising model in a transverse magnetic field, this KW duality symmetry 
maps the weak-field (low-temperature, in 2D context) ordered phase onto the strong-field 
(high-temperature) paramagnetic phase and vice versa. Although the strong-field phase 
appears disordered, it in fact sustains a hidden order which is revealed by a disorder operator 
0. Since the disorder operator is usually non-local and dual to the standard Ising order 
parameter, the two phases, which are separated by the Ising critical point and characterized 
respectively by the order- and disorder operators, are in many respects rather different from 
each other. The Ising (Z 2 ) quantum phase transition that occurs in this model can then be 
interpreted as a result of the competition between these two very different gapful orders Q. 

In this letter, we shall investigate several examples of one-dimensional competing orders 
whose critical properties are described, in the continuum limit, by a generalization of the 
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quantum sine-Gordon model with a manifest self-dual symmetry. The Hamiltonian density of 
the model is defined by: 

^SDSG = ^ 



'. + : cos (V8tt a r -Q^j '. J , (1) 

where the summation for r is taken over the positive roots of SU(AT) normalized to unity: 
6% = 1, and : : denotes the normal ordering symbol. The bo sonic vector field $ is made of 
iV — 1 free boson fields $ a ($ = . . . , $jv-i)) which are defined by chiral components 
$aR,L as: $ a = $ q l + $ a R, (a = 1, . . . , AT — 1). Similarly, each component of the dual vector- 
field 6 = (6i, . . . , 6jv_i) is defined by: a = $ aL — $ a R. The model (|T|) is a generalization 
of the usual sine-Gordon model where we have not only cosines of <3> but also those of the 
dual field 0. This field theory has been introduced in Ref. [4 J for exploring critical properties 
of vectorial Coulomb gas models in the presence of both electric- and magnetic charges. The 
interacting part of the model (UJ) is marginal and invariant under the Gaussian duality: $ <-» 
(i.e. the exchange of electric- and magnetic charges in the Coulomb gas context). In fact, 
as will be shown later, it has a hidden SO(iV) symmetry. Nevertheless, in what follows, the 
model © will be referred to as the SU(iV) self-dual sine-Gordon (SDSG) models. 

This model is of a direct relevance to the problem of competing quantum orders in 
one dimension. Indeed, since 0-field is a spatial integral of d t &, two fields $ and are 
mutually non-local and the model <[TJ) describes, in analogy with the above Ising duality, the 
competition between two completely different orders. In this respect, we shall give later 
several applications of the model (UJ)- For instance, one can anticipate that it describes the 
competition between a generalized charge-density wave, corresponding to the vertex operator 
of the $ field in Eq. (HJ), and a superconducting instability due to the perturbation depending 
on the dual field. The exact self-duality symmetry of the model <U]) may suggest the existence 
of a non-trivial quantum criticality in the infrared (IR) limit that results from this competition. 
In the simplest case (N = 2), the situation is well understood and a Gaussian U(l) criticality 
emerges whatever the sign of the coupling constant g J5l|6|. This model appears in the problem 
of the one-dimensional Fermi gas with backscattering and spin-non-conserving process as in 
the spin- 1/2 XYZ Heisenberg chain Q and also it describes critical properties of weakly- 
coupled Luttinger chains 0. The low-energy property for N > 2 is less clear. The 
perturbative study of the model O has been done in Refs. [SJIEl and a fixed point has been 
found whose nature has not been fully identified in these references. 

In this letter, we shall show that for g < 0, the model ([TJ displays a quantum critical 
behavior of the level-2 SO(iV) Wess-Zumino-Novikov-Witten (WZNW) universality class 
(hereafter the level-A; of the Kac-Moody algebra will be denoted as Gk) with central charge 
c = AT — 1. In contrast, for g > 0, it has a fully gapped spectrum and is related to an integrable 
field theory introduced by Andrei and Destri 0. 

The starting point of the solution is the introduction of A^ right-left moving Dirac 
fermions ^ a R,Li a = 1, . . . , N with free-Hamiltonian density: 

Ho = -i^ aR d x ^ aR + i^ aL d x ^ aL , (2) 
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where the summation over repeated indexes is assumed in the following. From these Dirac 
fermions, one can define SU(iV) "spin" chiral currents through: 

J R(L) = : ^I,R(L) T a/3 */3R(L) '■ , (3) 

where T A ,A = l,...,iV 2 — 1 are the generators of the Lie algebra of SU(iV) in the 
fundamental representation and normalized according to: Tr (T A T B ) = S AB /2. As well 
known, these currents satisfy the SU(iV)i Kac-Moody algebra and one can rewrite the free 
Hamiltonian © as a bilinear of currents (the so-called Sugawara form) iTTlfTOl: 



H = H 0c + H 0s = -(:Ji: + : Jt ) + { : J£J£ : + :J£J£\), (4) 

where we have introduced the U(l) "charge" currents: Jr(l) = ■ ^Li(L)^aR(L) : • At me 
level of the free theory Ho, spin and charge degrees of freedom decouple and the free "spin" 
Hamiltonian H 0s is nothing but that of the SU(iV)i WZNW conformal field theory (CFT). 
Note that the central charge of the model 7i 0s is c = N — 1, i.e. the central charge of N — 1 
massless free bosons which describes the g — > limit of the SDSG model (IB- 

Now let us add a perturbation V to the "spin" (or SU(iV)) Hamiltonian 7i 0s so that the 
"spin" part 7i 0s + V = Hn coincide with the sine-Gordon model dU). Obviously, it should 
be marginal (i.e. four fermion interaction) and invariant under both chiral (R<-s>L) symmetry 
and the Gaussian duality $ <-> 9. This self-duality symmetry considerably restricts the form 
of the possible four fermion interactions. To see this, let us introduce N chiral bosonic fields 
V^qR.l using the Abelian bosonization of Dirac fermions 0: 

* a R = —j== : exp (zV47r (p aR ) : (5) 

a / sir > ' 



/2tt 

^ah = —7= '■ exp (-iVATT (p aL ) : , (6) 
\J Itx v ' 

where the bosonic fields satisfy the commutation relation [(p a -R, = i$ a p/^- The 
anticommutation between fermions with different indexes is realized through the presence 
of Klein factors (here Majorana fermions) K a with the following anticommutation rule: 
{K a , Kp] = 25 a p. The Gaussian duality symmetry: ip a (= ip ah + (p aR ) <-> i? a (= <p aJj - ip aK ) 
thus amounts to the particle-hole (P-H) transformation only in the right-moving (R) sector of 
the Dirac theory: — > & aR , * aL —> ^ aL . As is well known, the SU(iV) generators T A 
can be classified into three categories: 

• Antisymmetric i.e. SO(iV) part: 

(^° (N) ) a/3 = ~(Wi)3 " (l<i<j<N) (7) 

• Symmetric part: 

(l*)^ = + (1 < i < j < N) (8) 

• Cartan generators (Diagonal): 

\ f m \ 

(r°) = = 51 *ofc*j9fc - "» <Wflfy,m+l I , (m = 1, . . . , A r — 1). (9) 

v 7a/3 ^2m(m+l) Vfc=i / 
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Since all generators belonging to the SO(iV) subset are antisymmetric, we deduce that 
the corresponding right currents behave under the Gaussian duality as: 

^-:»kW w 1.».*:=JES ,, \ 00) 



P R 1*3 /a/3 aR ■ ~ K ^ ' 

whereas the remaining (A r + 2)(A r — 1)/2 SU( N) generators are all symmetric or diagonal and 
the corresponding right currents change sign under the Gaussian duality: J R ' — > ~Jr ■ 
In contrast, the Gaussian duality does nothing for j£. This argument suggests that a possible 
model equivalent to the SDSG model <[TJ might be 

rtiV ~~ atTT ^ v • R R ■ + ■ L L ■ ) + ^ R ^' L>*i 

JV + 1 AeSU(N) i<j 

= £ £ [: (4T>) 2 : + : (^° (N) ) 2 : ] + A £ J L S °' N ', (1 1) 

for an appropriately chosen coupling constant A. In fact, by using bosonization rules Q-©, 
we can derive the SDSG model CO) from (fTTTl . Plugging Eqs. ©, © and © into (fTTTl . one 
obtains 



- A E { : cos(V5F(^ - ^)) : + : cos(v^F(7?i - ^)) : } . (12) 

If we introduce a charge bosonic field ^r^ and the SU(iV) bosonic fields $ r : l = 
1,...,JV- l)as El: 

$ C R(L) = Ol + • • • + ^Jv)r(l) 

$aR(L) = , (Vl + • • • + fa ~ ^a+lW) » ( 13 ) 

fa(a + l) 



the non-interacting part of Eq. (fT2l takes the standard form of a kinetic term for free bosons 
and the Hamiltonian (fT2l reads 



^ {: cos (v^ar-^) : + : cos (v^^r-©) : } -(14) 

r€A+ 



Thus we have shown that the model (fTTT) is indeed equivalent to the SDSG model O if we 
identify A = 87r 2 g. 

To deduce the physical properties of the model (fTTTl . it is more enlightening to introduce 
2N Majorana fermions £ RL and Xrl = !>•••) -^0 from the Dirac ones: Vl/jR^ 

30(N) 
l(L),ti' 

of the Majorana fermions: 



(£r,l + * Xr l) / V2. The ^r^w,-, being bilinears of Dirac fermions, can be expressed in terms 



J R(L),jj _ 2 V^R( L )^R(L) ^ A-R(L)AR(L)J — 2 R ( L )«' 
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where J^AJj, being the sum of two SO(iV)i currents, is an SO(iV)2 current. Therefore, we 
deduce that the SDSG model © is equivalent to the level-2 SO(iV) WZNW model perturbed 
by a marginal current-current interaction: 

Wsdsg = £ £ [ : PSD' : + : (JiT? ] + ^ £ tfTtfT- (W 

This equation is one of the main results of this letter. 

Using this equivalence, one can extract the IR properties of the SDSG model. The one- 
loop renormalization-group (RG) equation of the model (fToT) is: g = (N — 2)Tig 2 , where 
g = dg/dl (I being the RG parameter). For g < 0, the interaction is marginally irrelevant so 
that in the far IR limit, the model flows towards the SO(iV) 2 WZNW fixed point with a central 
charge c = N — 1. This CFT has the same central charge as the SU^AOi WZNW model and in 
fact there exists a conformal embedding between them ifTUll: SU(A^)! D SO(iV) 2 . In contrast, 
when g > 0, the interaction is marginally relevant and flows toward strong coupling. From 
the structure of the current-current interaction, it is naturally expected that a mass gap opens 
dynamically i.e. the SDSG model is a massive field theory for g > and N > 2. In fact, this 
can be explicitly shown by observing that model (fToT) is related to an integrable field theory 
introduced by Andrei and Destri [9| (see also Ref. C3) with the following Hamiltonian: 

Had = -~ l x dx^\,i + $2,i ^dx^i) - 9ad(p 2 + p 2 + a 2 + cr 2 ), (17) 

where i^n (respectively i/> 2 ,i) is a two-component spinor formed by ^ L (respectively Xr l) 
and 7° = a 2 , 7 1 = io\, and 75 = a 3 (<Tj being the Pauli matrices). The 0(iV)-invariant order 
parameters p, a, p and a are defined as: 

p = \ + = -i + xkxi) (is) 

a = - farfikv = -i (^rXL - Xitl) (19) 

P = \ - M2,i) = -1 - xixi) (20) 

5 = ^2,i = -i(^X*L + XR^)- (2D 
Although Had looks complicated, after bosonizing, the Hamiltonian <TTvT> separates into two 
commuting pieces, a free Hamiltonian 7i 0c for the massless bosonic field $ c (Eq. (TH51) ) and 
the SO(A^) 2 current-current model (fT6b with g = ^ad/tt 2 : Had = Hoc + ^sdsg- Since 
the model ( fT71) is exactly solvable by means of the Bethe ansatz we can extract the 

physical properties of 7Ysdsg from the solution. The nature of the ground states may be simply 
understood in terms of the order parameters p, a, p and a, which form two independent SO(2) 
doublets (p, a) and (p, a) These two doublets are mapped onto each other by the 
KW duality for the Majorana fermions x u - Xr ~^ ~Xr an d Xl ~^ Xl- F rom the form of the 
interacting part of the model (fTTt . we readily see that it is invariant under the interchange of 
the two doublets. On the basis of large- N semiclassical argument, the authors of Ref. |j9| 
found that when #ad < 0, this interchange symmetry is broken spontaneously in the ground 
state and that there are two different ground states where only one of the two doublets has a 
finite modulus; correspondingly massive kink excitations appear in the spectrum to connect 
the above ground states. 
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Let us now consider some physical applications of the SDSG model © in the context of 
competing orders in one dimension. As the first example, we take a spin-1 bilinear-biquadratic 
Heisenberg chain lfT4ll with nearest ( J\) and next-nearest ( J 2 ) interactions: 

2 

Hbb = ^ ^ Jbb S„ • S n+a + (S n • S n + a ) + 5bb ^ S n • S n +i, (22) 

n a=l n 

with S n being a spin-1 operator at site n. The model with 5bb = is SU(3)-symmetric 
and, in particular, for Jgg = it reduces to an integrable model lfT31 fT6l which displays a 
quantum critical behavior of the SU(3)i WZNW universality class ifTTl [T8l . The effect of 
the remaining interactions can be investigated in the vicinity of the SU(3) symmetric point 
= 5bb = 0) using the low-energy approach of Ref. llT8l . In fact, Itoi and Kato 
|[T8l considered a more general problem of an SU(iV)i WZNW model perturbed by the most 
general SO(iV)-symmetric marginal perturbation: 

Wnc = j^j ( : 44 ■ + : 44 : ) + ^44 + 2A 2 (t*T*) (23) 

which for N = 3 should describe the low-energy physics of the SO(3) model d22t around the 
SU(3) symmetric point. 

Using the decomposition (171819b of the SU(iV) generators T A , the model (l23l) can be 
expressed in the following compact form: 

Wnc = t^t + : 44 ■ )+(Ai - A 2 ) £ + A 2 ) J2 44 m) 

iV + 1 agso(n) Aes,D 

It is straightforward to calculate the one-loop RG equations lfT8ll for the model (l24b and we 
obtain 

G\ = — G\ H — G\ , G 2 = —GiG 2 , (25) 

671 8tt 4tt 

where we have introduced a new set of couplings as G\ = \± — A 2 and G 2 = Ai + A 2 . The 

RG-flow is shown in Figure [T] We have two gapful phases (phase-1 and phase-2) together 

with one extended gapless phase which belongs to the SU(iV)i WZNW universality class. 

In the special case of N = 3, the Hamiltonian d24b describes the competition between 

two gapful orders of i^BB: a trimerization (period-3) phase (phase-1 in Figure [T]) stabilized 

when A 2 = and Ai > 0, where three adjacent spins form local SU(3) singlets, and the 

non-degenerate Haldane state (phase-2) when Ai = and A 2 < lfl8lfT9l . The trimerized 



(2) 

phase is expected to occur in Hbb when 5bb = and for a sufficiently strong value of J BB \ 
whereas the Haldane phase appears when J B g = and 5bb > 1141 . The one-loop RG 
flow is presented in Figure [T] and we see that the phase transition between these two gapful 
phases occurs along the line Ai = — A 2 > shown as 'SDSG'. From Eq. d24l . we find that 
the effective field theory which describes the transition is given by the SO(iV) current-current 
model <n~6t with N = 3 and g = \i/4it 2 (= — A 2 /47r 2 ) > 0. We thus have found an example 
of the SDSG model ^ with N = 3 which describes the competition between the Haldane and 
trimerized orders and corresponds to a first-order transition (since the gap opens for g > 0). 
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Figure 1. One-loop RG flow for the Itoi-Kato model (1231 . In the sine-Gordon language, 
the Ai- and the A2 axis respectively correspond to the pure cos(y/8nd r -<^) model and the 
pure cos(v / 87ra?r @) one (see Eq. Q). Since they are related to each other by the Gaussian 
duality, phase-1 and phase-2 cannot be connected by any local symmetry. These competing 
gapful phases are separated by the line of self-dual sine-Gordon model (SDSG). The thick 

(2) C2l 

arrow schematically shows the path traced by H^(J^ = 0) (H S0 (Js = 0) for N = 4) as 
Sbb (<5so for N = 4) is changed from positive (Haldane- or staggered dimerization phase) to 
negative (gapless SU(./V)i WZNW phase). 

In the second example, we consider an SU(2)xSU(2)-symmetric spin-orbital chain with 
nearest (J^) and next-nearest ( J@>) interactions: 

#so = EE J so ] (2 S n -S n+a + I-) (2 T n -T n+a + 5 so £ (S n -S n+1 + T n -T n+1 ) , (26) 

where S„ and T„ denote spin- 1/2 operators representing respectively the spin- and the two- 
fold degenerate orbital degrees of freedom ll2"Ul ITU on the n-th site. For jj£> = S so = 0, 
the model coincides with an SU(4) generalization of the spin- 1/2 Heisenberg chain and is 
exactly solvable by Bethe ansatz [ 16] ; the model is gapless with three massless bosonic modes 
and the field theory describing this quantum criticality is the SU(4)i WZNW model lfT71l22l 
or, equivalently, the SO(6)i WZNW theory in terms of two triplets of Majorana fermions 
CrliXrli* = 1, • • • , 3 ll23l [T3ll . Using the Majorana basis, one can derive the low-energy 
effective Hamiltonian of the model d26T) in the vicinity of the SU(4) point (Jj$ = 5 so = 0). 
We find, using the results of Refs. lfT3ll23l . the following effective Hamiltonian density: 

+ (9i + 92) (CkCl + XrXI) 2 + (91 - 92) (CkCl - XkXl) 2 ■ (27) 
Since the interacting part of Eq. d27b can be written as: TC^ = —{91 + 92) P 2 — (pi — 5 , 2)p 2 , 
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we may expect that the model describes the competition between two different fully gapped 
orders which are characterized, in this continuum limit, by the order parameters p and p of 
Eqs. (II 81201) . As has been shown in Ref. lfT3ll . they correspond respectively to an SU(4) 
quadrumerization (Q) phase (i.e. period 4), which is characterized by local SU(4) singlets, 
and to a period-2 staggered dimerization (SD) phase which is formed by alternating spin and 
orbital singlets ll24l . In terms of the lattice coupling constants, the former phase emerges when 
5 so = and for a sufficiently large value of Jj£' ll2"51 whereas the SD phase is stabilized when 
J$ = and 5 so , jW > EDI . The competition between these two orders can be investigated 
by observing that the SU(2)xSU(2)-symmetric model H so (l26l) can be recasted into the form 
of Hik (EH) with N = 4 since SU(2)xSU(2)~SO(4). In fact, an explicit calculation shows 

= §51 E J R J l + % 2 E J rJl ■ (28) 

AeSO(4) AeS,D 

With the identification G\$ = 8g 1>2 and N = 4, the RG equations d25t again describe the 
model H so . From the Figured we observe that the quantum phase transition between the two 
competing phases (phase-1 for 'Q' and phase-2 for 'SD') occurs at g 2 = and g± > and is 
thus described by: 

H so {g 2 = 0) = 7Yik(Ax = -A 2 = A 9l ) = H^ G (g = 9l /n 2 ) . (29) 

From the equivalence between 7Ysdsg an d the massive sector of Had, we conclude that the 
Q<-*-SD transition described by H so (g2 = 0) is of first order. 

The last example is provided by the generalized two-leg spin ladders with four-spin 
exchange interactions studied recently in Ref. lfT31 . In particular, it has been shown that, close 
to the SU(4) symmetric point of Eq. d26T) with jj£> = S so = 0, four competing orders emerge. 
In addition to the Q (p) and SD (p) orders of the previous example, a scalar-chirality order lf26l 
and a rung-quadrumerization order appear. These two additional phases are characterized, 
within the low-energy approach of Ref. lfT31 . respectively by the order parameters a and a of 
Eqs. (I19I21I) . As is seen from the one-loop RG analysis of Ref. fT31 . the competition between 
these four orders is governed by the low-energy effective Hamiltonian: 

^cfr = -y (&d x & - &d x & + XrQ^r ~ Xl^xl) - A (p 2 + p 2 + a 2 + a 2 ) , (30) 

with A > 0. We thus observe that the Andrei-Destri model (fTTt with N = 3 and g^D = A 
accounts for the competition between the four different orders of the problem. Since the latter 
model is equivalent to the SDSG model (OQ) up to a free massless bosonic field, we easily see 
that the resulting phase transition is of a U(l) Gaussian type when A > 0. 

We hope that other applications of the SDSG model will be reported in the near future. 
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